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The invention relates to symbol lock detector for coherent digital 
communication systems after carrier and subcarrier, if any, synchronization 
has been achieved in a phase-locked loop employing nonoverlapping symbol 
intervals, where symbol data, d^, takes on binary values with equal 
probability, and symbol interval T is known. 

FIG. 1 is a generic functional block diagram, while FIGs. 2a and 2b 
are functional block diagrams of respective prior-art square-law symbol 
lock detector (SQOD) and prior— art absolute value symbol lock detector 
(AVOD), with integration of overlapping intervals in both cases, thus 
requiring separate integrators 21 and 22 in both cases. FIGs. 3a and 3b 
are functional block diagrams of a new nonoverlapping square-law symbol 
detector SQNOD and a new nonoverlapping absolute value detector AVNOD, 
respectively, that process nonoverlapping half symbol intervals (T/2) so 
that a single integrator 30 may be used in each case with a multiplexer 32 
for routing the inphase data through a half symbol interval delay element 
33 to means for forming a difference X k = I k 2 -Q k 2 and accumulating a number, 
M, of X k samples to form a statistic Y for threshold detection, all with 
functional blocks 26, 27 and 28 in series. A third new signal— power symbol 
lock detector (SPED) shown in FIG. 4 omits both the squaring operation 31 
and the absolute value operation 31' of the embodiments of FIGs. 3a and 3b 
and instead utilizes a multiplexer 41 to connect the inphase data I k 
through a delay element 42 to a multiplier 43 and to connect the quadra- 
phase data Q k directly to the multiplier 43 to form a product X k over many 
symbol intervals and thresholding an accumulation Y of AT samples of the 
product X k to form a symbol lock decision when the threshold 6 is exceeded. 
This provides greater simplicity in the implementation of a symbol lock 
detector with better performance when the threshold <5 is set in the pre- 
sence of noise only (no signal) than even the prior-art symbol lock detec- 
tors which outperform the first two new symbol lock detectors of FIGs. 3a 
and 3b when the threshold is set in the presence of a signal. The first 
two new symbol lock detectors nevertheless provided sufficiently good per- 
formance for them to be considered for DSN communication receivers because 
of their advantage of requiring only one integrator. FIGs. 5a and 5b 
through FIG. 8 present graphs of analytical data which uphold these 
conclusions about performance. 

The novelty of the invention resides in reaching a symbol lock deci- 
sion by integrating nonoverlapping half symbol intervals so that only one 
time-shared integrator is required. An electronic multiplexing switch re- 
quired in order to operate with just one integrator is very simple as com- 
pared to implementing a second integrator. A further improvement resides 
in forming the product (power) of the integrals of nonoverlapping half sym- 
bol data instead of first forming the squares or absolute values and then 
forming the sum of their differences. 
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SYMBOL LOCK DETECTION IMPLEMENTED 

WITH NONOVERLAPPING INTEGRATION INTERVALS 

ORIGIN OF INVENTION 

The invention described herein was made in the 
performance of work under a NASA contract, and is 
subject to the provisions of Public Law 96-517 (35 USC 
5 202) in which the contractor has elected not to retain 

title . 

TECHNICAL FIELD 

This invention relates to symbol lock detectors 
for coherent digital communication systems after 
10 carrier and subcarrier synchronization has been 
achieved in a phase lock loop, and more particularly 
to lock detectors which may be implemented with a 
single integrator for processing nonoverlapping symbol 
intervals of an incoming signal. 

15 BACKGROUND ART 

Deep space network (DSN) receivers currently 
under development use phase locked loops (PLLs) to 
track the carrier, subcarrier, and symbol phase. Like 
most coherent receivers, the DSN receivers rely on 
20 lock detectors to provide the symbol lock status of 
its PLLs. Since carrier, subcarrier, and symbol 
synchronization need to be achieved before any mean- 
ingful symbol detection can be initiated, symbol lock 
detectors play a vital role in the final decision of 
25 accepting or rejecting the detected symbols. In the 
past, symbol lock detectors have employed overlapping 
symbol intervals in their operations and therefore 
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require two integrators operating over staggered time 
intervals. 

During operation, a loop is assumed to be locked 
when its lock indicator consistently has a positive 
5 status. The carrier and subcarrier lock detectors 
currently used in DSN receivers have already been 
analyzed. The present invention concerns the analysis 
of three new symbol lock detectors which simplify 
implementation for DSN receivers and two prior-art 
10 symbol lock detectors for comparison. 


Statement of th e Invention 

In accordance with the present invention, an 
incoming signal that has undergone carrier and subcar- 
rier (if any) synchronization is processed to deter- 
15 mine whether or not symbol lock has also been 

achieved. The symbol data, d k , takes on the binary 
values of a pulse p (t) of a known duration, such as in 
a nonreturn-to-zero (NRZ) or Manchester code, where 
the probability of the data having a +1 is equal to 
20 having a -1 (or 0) value, (i.e., the data transition 
probability equals one half). While the prior art 
lock detectors employ two integrators for processing 
two overlapping symbol intervals, the lock detectors 
of the present invention require only one integrator 
25 for processing two nonoverlapping time intervals, 
e.g., integration of the last half of one (<W 2) and 
the first half (d*/2) of the next symbol interval, 
followed by integration of the second half (d k / 2) of 
the next and the first half (d k+1 /2) of the following 
30 symbol interval, or the first half of each symbol 
interval followed by the second half of each symbol 
interval. The receiver is assumed to have perfect 
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knowledge of the symbol interval T=kt, where k is the 
number of data bits, and t is the data bit interval. 

The first integration output is processed and de 
layed a half symbol interval, T/2, for determining 
5 symbol lock by addition to or multiplication with the 
following integration output, where the processing of 
the first and second integration outputs is a squaring 
operation or an absolute value operation, and symbol 
lock is determined by addition of the processed second 
10 integration output to the processed first integration 
output delayed a half symbol interval and accumulating 
the sums. Symbol lock is then determined by repeated- 
ly thresholding the accumulated sums every M symbol 
intervals. An alternative arrangement for determining 
15 symbol lock is by multiplying the second integration 
outputs by the first integration output delayed a half 
symbol interval to obtain a signal power estimation, 
and after accumulating signal power estimates over M 
symbol intervals, thresholding the accumulated esti- 
20 mates. 

BRIM, description of the drawings 

FIG. 1 is a functinal block diagram common to 
symbol lock detectors shown in FIGs. 2a through 4. 

FIG. 2a is a functional block diagram of a prior- 
25 art square-law symbol lock detector with integration 
of overlapping intervals (SQOD). 

FIG. 2b is a functional block diagram of a prior- 
art absolute-value symbol lock detector with integra- 
tion of overlapping intervals (AVOD) . 

FIG. 3a is a functional block diagram of a 
square-law symbol lock detector with integration of 


30 
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nonoverlapping intervals (SQNOD) in accordance with 
the present invention. 

FIG. 3b is a functional block diagram of an 
absolute-value symbol lock detector with integration 
5 of nonoverlapping intervals (AVNOD) in accordance with 
the present invention. 

FIG. 4 is a functional block diagram of a signal- 
power symbol lock estimator (SPED) with integration of 
nonoverlapping (first and second) halves of incoming 
10 signal symbols in accordance with the presnet inven- 
tion . 

FIG. 5a and 5b are graphs of the probability 
density function of the output Y for the SQOD of FIG. 
2a when symbol lock is not present and the output Y 
15 has a high SNR=5dB and a low SNR=-5dB, respectively. 

FIG. 6 is a graph of the probability of lock 
detection when r is an unknown constant over a deci- 
sion interval for symbol lock detection of FIGs. 2a , 
2b, 3a, 3b and 4. 

20 FIG. 7 is a graph of the probability of lock 

detection versus SNR when r is uniformly distributed 
and changing from symbol to symbol for symbol lock 
detectors of FIGs. 2a, 2b, 3a, 3b and 4. 

FIG. 8 is a graph of the probability of lock 
25 detection versus SNR when the false alarm rate is 
computed in the absence of a signal for symbol lock 
detectors of FIGs. 2a, 2b, 3a, 3b and 4. 

DETAILED DESCRIPTION OF THE INVENTION 

The symbol lock detectors considered are divided 
30 into two groups. The detectors in the first group are 
prior-art symbol lock detectors that process the 
overlapping outputs of two symbol data bit integra- 
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tors, whereas those in the second group, which depart 
from the prior art, use one integrator for two non- 
overlapping intervals of an incoming signal to deter- 
mine whether or not symbol lock has been achieved. 

5 The first group of symbol lock detectors employ 

either the squares of the two overlapping integrator 
outputs or the absolute value of the integrator 
outputs. The second group consists of three new lock 
detectors, two of which use the same mathematical 
10 operations of squaring or taking the absolute values 
of integration outputs of nonoverlapping half symbol 
intervals but only one integrator, while the third de- 
tector, which also uses only one integrator, does not 
include either squaring or absolute value operations, 
15 and instead functions as a signal power estimator by 
multiplying the two integration outputs of nonoverlap- 
ping half symbol intervals. 

The five lock detectors are compared based on the 
lock-detection probability as a function of the symbol 
20 SNR for a given false-alarm probability and a fixed 
observation interval. 

Although symbol lock detection has been addressed 
before [J.K. Holmes, Coherent Spread Spectrum System, 
New York: John Wiley and Sons, 1982 and K.T. Woo, 

25 Shuttle Bit Synch Lock Detector Performance, TRW IOC 
No. SCTE 50-76-184/ KTW , TRW Corporation, El Segundo, 
California, April 5, 1976], the analyses have ne- 

glected the interdependence between symbol synchro- 
nizer bandwidth, and lock detector bandwidth. The 
30 symbol synchronizer bandwith refers to the one-sided 
loop noise bandwidth B L of the digital data-transition 
tracking loop [M. Simon, "An Analysis of the Steady- 
State Phase Noise Performance of a Digital Data- 
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Transition Tracking Loop," JPL Space Programs Summary, 
37-55, Vol. 3, Jet Propulsion Laboratory, Pasadena, 
California, pp. 54-62, February 28, 1969] used in the 
receivers. The lock detector bandwidth is defined as 
5 the frequency at which the lock detector provides a 
status signal of being in or out-of-lock. For exam- 
ple, the lock detectors considered here indicate loop 
status once every M symbols. Consequently, the 
bandwidth of these detectors is 1/MT, where T is the 
10 symbol interval. 

The probability of false alarm, Pj a , is defined in 
two ways. In the classical sense, it is defined as 
the probability of declaring a signal (or target as in 
radar applications) to be present when it is not 
15 present. In deep space applications however, it is 
more appropriate to define Pj a as the probability of 
declaring a loop to be in-lock when it is out-of-lock. 
That is, declaring the timing error to be zero (in- 
lock) when the loop is slipping cycles and operating 
20 with a non-zero timing error (out-of-lock) . 

In discussion below, the false alarm rate is 
shown to be drastically different depending on the 
definition used. In addition when the loop is slip- 
ping cycles, the false alarm rate is shown to depend 
25 strongly on the ratio of the lock detector bandwidth 
to the symbol loop bandwidth. For example, when the 
loop is slipping and 1 /B L « MT the lock detectors 
operate with acceptable false alarm rates because 
there are several uncorrelated samples of the timing 
30 error 1/r within the MT second decision interval. On 
the other hand, when 1 /B L » MT the false alarm rates 
are unacceptable because the timing error is constant 
over several decision intervals. Note that a good rule 
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of thumb is to assume that the loop provides uncorre— 
lated phase estimates every 1/B l seconds. As a 
result, the symbol timing error at time t i is uncor- 
related with the symbol timing error at time tj when 
5 I fi “ tj I > 1 /B l . This article considers the special 
cases of 1 /B L =MT and 1/B L = T. The first case is ana- 
lyzed and simulated whereas, the second is simulated 
but not analyzed. When the threshold is adjusted in 
the presence of noise only, the performance can be 
10 derived from the previous analysis by setting the 
signal amplitude to zero. 

Generic Description of Lock Detectors 

FIG. 1 is a block diagram showing the signal pro- 
cessing functions common to the symbol lock detectors 
15 analyzed and discussed below, including both prior-art 
lock detectors of FIGs. 2a and 2b, and new lock detec- 
tors in accordance with the present invention as shown 
in FIGs. 3a, 3b, and 4. The received signal is assumed 
to have been mixed with perfect carrier and subcarrier 
20 local reference signals so that the input to the lock 
detectors is a baseband signal of the form 

r(t) = Adit) + n(t) (1) 

where A is the signal amplitude and A 2 is the received 
data power with 

d(t) = dy-pit-kT) 


( 2 ) 
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and where n(t) is the additive white Gaussian noise 
process with a single-sided power spectral density 
(PSD) N 0 (Watts/Hz). The data symbol d k takes on the 
binary value +1 and -1 (or 0) with equal probability, 

5 and p(t) is the received data pulse shape of duration 
T seconds. For comparison purposes, only the NRZ pulse 
is considered in the analysis but the results can be 
extended to any pulse shape, such as Manchester 
encoded data. The receiver is assumed to have perfect 
10 knowledge of T, but not the symbol epoch, i.e., the 
receiver has estimated perfectly the symbol rate but 
not necessarily the start and end of the symbols. 

The signal processing functions for the symbol 
lock detector 10 in FIG. 1 depend on the processing of 
15 integration outputs in the detector. Its output X k is 
at the symbol rate and typically many samples of X k 
are averaged in an accumulator 12 to obtain the 
decision statistic Y compared in a block 14 with a 
threshold value 6. If Y is greater than the threshold 
20 6, the loop is declared to be in— lock, otherwise it is 

declared to be out-of-lock. 

The symbol timing error (parameter) r in FIG. 1 
is the phase error between the symbol phase and the 
phase estimate provided by the symbol synchronizer. 
25 The in-lock case is analyzed by setting the timing 
error t to zero. In practice, the error is not 
identically zero, but it is a very small value. When 
there is a signal present, the out-of-lock model for 
t depends on the relation between B L and 1 / (MT) . When 
30 B l = 1 / (MT) , r is modeled as an unknown constant over 
a decision interval (MT seconds) but independent and 
uniformly distributed from one decision interval to 


the next. 
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Alternatively, when B L = IT, the timing error is 
modeled as constant over a symbol interval, T, but 
independent and uniformly distributed from symbol to 
symbol. In this case, if the decision time MT » T (as 
5 xt usually is), each decision statistic encompasses 
the entire range of r. When there is no signal pres- 
ent, the model of r is irrelevant because the out-of- 
lock performance is independent of r. 

The respective prior-art and new detectors consi- 
10 dered and compared below are the square-law detector 
with overlapping (SQOD) and non-overlapping (SQNOD) 
integrators shown in FIGs. 2a and 3a, the absolute- 
value detectors with overlapping (AVOD) and non-over- 
lapping (AVNOD) integrators shown in FIGs. 2b and 3b, 
15 and finally the new signal power estimator detector 
(SPED) shown in FIG. 4. 

In the prior-art SQOD detector, FIG. 2a, the in- 
put signal r(t; is integrated over two symbol periods 
by two integrators 21 and 22: one in phase with the 
20 estimated symbol interval and the other staggered by 
half a symbol duration. The resulting inphase and 
quadraphase samples I k and Q k are correlated due to the 
overlapping intervals. The quadrature samples I k are 
squared in a processor 23 and delayed a half symbol 
25 period by a delay element 24 while the samples Q k are 
squared in a processor 25. The delayed ( l k ) 2 and the 
undelayed (Q k ) 2 samples are then combined in a summing 
circuit 26 to form an output X k which are averaged in 
an accumulator 27 to obtain the decision statistic Y 
30 which is then compared with a predetermined threshold 
6 in a detector 28 to reach a lock decision. The 
prior-art AVOD detector, FIG. 2b merely replaces the 
squaring operations in FIG. 2a with absolute value 
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operations in blocks 23' and 25'. Hence, the I and Q 
samples of the prior-art AVOD detector are also cor- 
related. 

The new SQNOD detector, FIG. 3a, processes the 
5 integration outputs of nonoverlapping half symbol 
intervals from a single integrator 30. A single 
squaring operation in block 31 produces the square of 
the integration outputs from both halves of a symbol 
interval, and a multiplexer separates the first and 
10 second half symbol interval integration outputs 
sending the first half through a delay element of half 
a symbol interval and the second half directly to a 
summing circuit 34. Its output X k is then averaged in 
an accumulator 35 and its output, and the decision 
15 statistic Y is compared with a threshold 6 in a de- 
tector 36. As before, replacing the squaring opera- 
tions in FIG. 3a with absolute value operations in 
block 31' yields its counterpart the AVNOD detector 
shown in FIG. 3b. The integrator outputs in these 
20 cases are uncorrelated because the integrated inter- 
vals are nonoverlapping. 

The SPED detector shown in FIG. 4, which also 
uses a single integrator 40, was considered for symbol 
lock detection because it already existed as part of 
25 a split symbol moment SNR estimator [K.T. Woo, supra] 
in DSN receivers used for SNR estimation. The inphase 
( I k ) and quadraphase (Q k ) integration outputs of the 
SPED are obtained by integrating the received signal 
over the respective first and second halves of a 
30 symbol interval and separating them by a multiplexer 
41. Since the noise in the first and second half are 
independent, delaying the inphase integration output, 
I k , in a delay element 42 and then forming the product 
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of I* and Q k in a multiplexer 43 and averaging over 
number, M, of symbol intervals in an accumulator 44 , 
provides as a decision statistic Y an estimate propor- 
tional to signal power which is then compared to the 
value 6 in a threshold detector 45 to provide a lock 

decision * 


Performance Analysis 

In-lock performance is measured in terms of the 
probability of declaring a phase-locked loop (PLL) as 
10 being symbol locked when there is no timing error, 
i.e., the probability that the decision statistic Y is 
greater than the threshold value 6 when z = 0. Note 
that r = 0, or no phase tracking error, is equivalen 
to setting the symbol synchronizer loop SNR to infim- 
15 ty The degradation in detection probability due to 
timing jitter (non-infinite loop SNR) is minimal and 
has been addressed in the case of carrier lock detec- 
tors [A. Mileant and S. Hinedi, "Costas Loop Lock 
Detection in the Advanced Receiver," TDA Progress 
20 Report 42-99 , Vol. July-September 1989, Jet Propulsion 
Laboratory, Pasadena, California PP . 72-89, November 
15, 1989]. The out-of-lock performance is measured by 
thl probability of false alarm, i.e., the probability 
of declaring the loop as locked when it is not locked. 

25 The out-of-lock performance in the presence of a 
signal is analyzed for the case B h = 1/(MT). Note 
that, in that case, the timing error r is independent 
from one symbol to another and the decision is per- 
formed after averaging a number, M, of symbols. On 
30 the other hand when = U (MT) , z is an unknown 

constant during a decision interval and independent 
from one decision to the next. Setting M • 1 xn the 
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latter case would imply a decision every symbol, which 
is fundamentally different from the case B L = 1/T 

decision using M symbols. Hence, the performance when 
B l = 1/T cannot be derived from the case of B L = 

5 1/ (MT) , simply by setting M = 1. The out-of-lock 

performance when there is no signal present is also 
analyzed. 

In this discussion, only the final equations are 
shown. Derivations of these equations are set forth 
10 in the various appendices A through F for all five 
symbol lock detectors SQOD, SQNOD, AVOD, AVNOD and 
SPED which, by this reference, are hereby made a part 
hereof. In all cases, the decision statistic can be 
expressed as the average of samples X k over M symbol 
15 intervals given by the equation 


Y = 


M Jt-1 


(3) 


Note that the random variable X k is peculiar to each 
detector. When the timing offset r = 0, the adjacent 
samples X k and X k+1 are correlated in the two prior-art 
detectors (SQOD and SQNOD), whereas for the three new 
20 detectors (AVOD, AVNOD and SPED) they are uncorre- 
lated. In all cases, the random variable X k is not 
Gaussian due the nonlinear operations on I k and Q k . 
For large M, the random variable Y is modeled as 
Gaussian due to the Central Limit Theorem (CLT) . The 
25 CLT theorem applies to the sum of correlated random 
variables when none of the variables being summed or 
multiplied dominates over the others [D. Fraser, Won- 
Parametric Methods in Statistics , New York: John 

Wiley and Sons, 1957]. This model for Y is justified 
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by simulation results. The probability of lock detec- 
tion is the probability that the Gaussian random 
variable Y surpasses the threshold 6. Hence, it is 
given by 


Pd 



exp 


(y-M 2 

20y 


dy 


( 4 ) 


5 where yy and Oy are the mean and variance of Y when r 

is exactly zero. Using the definition of the error 
function 


erf(x)k~ f* exp ( -t 2 )dt 


one has 


Pd 


1 

2 


1 

2 


erf 




( 6 ) 


or 


1-A erf 

6 

»' " \ 

snr d 

2 2 

^20y > 

2 

/ 


(7) 


10 where SNR D denotes the detector SNR defined by 

srn, « (fj 


( 8 ) 
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The threshold 6 is chosen to maintain a fixed 
probability of false alarm. The probability of false 
alarm is the probability that the out-of-lock decision 
statistics do not surpass the threshold. Hence, it is 
5 given by 

p f a = f 0 (y) d y <9) 


where / Q (y) is the out-of-lock density of Y . The 
threshold 6 is computed by solving Eq. (9) for a fixed 
Pj a . When there is a signal present and B L = 1/(MT), 
the statistic Y is no longer Gaussian and f Q (y) must 
10 be obtained numerically or by simulation as shown 
below. When there is no signal present, the CLT can 
be invoked and the out-of-lock decision statistic can 
be modeled as Gaussian. This model is verified by 
simulations below. In this case, Eq. (9) can be 
15 written as 


/a 


1 

2 


1_ 

2 


erf 


flltai 

\fZ°r 0 


( 10 ) 


where yi yo and a yo the out-of-lock mean and variance of 
the decision statistic Y . The threshold 6 is given by 


6 = v/2o yo v + H ro 


( 11 ) 


where y = erf _I (l - 2 P fa ) • Substituting Eq. (11) into 
Eq. (6) relates the probability of detection to the no 
signal (classical) false alarm rate, namely. 


20 
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I _ A erf f v ^ qv ° Y + >ly ° — 
22 ^ /2o r 


( 12 ) 


The next five subsections derive the in-lock and 
out-of-lock mean and variance for all five detectors 
SQOD , SQNOD, AVOD, AVNOD and SPED for comparison. 

(1). Square-Law Lock Detector with Overlapping Inter 
5 vals (SQOD) 

The SQOD detector is shown in FIG. 2a. For the 
input given by Eq. (1), the inphase integrator output 
is given by 


I* = 


f Uc'l 

J Jtr** 


(ku) r*t 


r (t) d£ 


(13) 


= d*A (T-x) +d ktl Ax+N l (k) +N 2 (k) 


and the quadraphase integrator output is given by 


/*(*+})T + r 

Qic — / r(t)dt 

J(H-i)T+r 


' dkA (y — r) 4- dk+iA (y + r) 

+w 2 (jt) + Wi(Jt + i) o<r<r 

+ (^* “ r ) + dk+iA (r — y) 

, + N 2 (k) + N,(k + 1) 


£ <t<T 
(M) 
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where r is limited to the interval [O f T] , and where 

( 15 ) 


N,(k) = |* 2 ' n(t)dt 

J kT+* 


and 


N 2 (k) 


-f, 


(Jc+l) T+ T 

( jc 4) T+t 


n{c) dt 


( 16 ) 


Since nft) is a white Gaussian process with one-sided 
5 PSD N ol the W^'s are independent Gaussian random vari- 
ables with mean zero and variance o^=(N 0 T)/ 4. The 

samples X k =ll~Ql and summing M of them yields Y. From 

Appendix A, the in-lock mean and variance of Y are 
given by 

„ _ ( 17 ) 


10 and 



( 18 ) 


The out-of-lock mean and variance when there is a 
signal present and r is an unknown constant over a 
decision interval {B L = 1 / (MT) ) are given by 
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Vy a = °-° 


(19) 


'Yo 


[ n 0 t\ 

/ 


M 

V M J 

. \ 


31r\ z e+ 41i] t 


120 


24 


+ -± +2(Af-l) 


23ri| 23 ti s 
320 48 



+ (M-1) ( M-2 ) 



( 20 ) 


where ri s denotes symbol signal-to-noise ratio and is 
defined as 


T\ a ± £1 

No 


( 21 ) 


5 


Setting A = 0 in Eq. (20) and substituting the result 
in Eqs . (19) and (20) yields the out-of-lock mean and 
variance in the no signal case. Hence, the no signal 
mean is zero but the no signal variance is given by 




¥f ["i- 2 *"- 1 ’-!] 


( 22 ) 


(2). Absolute-Value Lock Detector with Overlapping 
Intervals (AVOD) 

For the AVOD detector shown in FIG. 2b, absolute 
values are used instead of squares. The expressions 


10 
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for I k and Q k given by Eqs. (13) and (14) are still 
valid but now X k = | I k | - 1(3*1 . From Appendix B, the 
in-lock mean and variance are given by 


= 


N a T 




exp (— n s ) -1 


( 23 ) 


erf (/id 


and 


a l = — [M Var(X k ) Cov(X k ,X k ^) ] ( 24 ) 

M 2 


5 where 


Var (X k ) 


(N q T) 



+l-i) s [F 1 (i) s ) +F 2 ( t\ s ) 


] 


[exp (-2r) s ) -2exp(-t) s ) +1] — ^ erf 2 (/r^ 
4tc 4 


-^erf (JrT s ) [exp(-ii s ) -1] J 


( 25 ) 


and 


Cov(X k ,X k tl ) = (N C T) 


exp (~2ti s ) +exp (-r) s ) 


2n 


> 
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erf 2 (y/tQ 


exp(-ti s ) 




[F^n,) *F 2 (tj s ) ] 


( 26 ) 


The out-of-lock mean and variance when there is 
a signal and tO is constant over M symbols are given 
by 


Vyo = °-° 


(27) 


and 


0 y o = -±-lMVar 0 (X k )+2(M-l)Cov 0 (X k ,X k . 1 ) 
M 2 


+ {M- 1) ( M-2 ) Cov 0 (X k , X k . 2 ) ] 


(28) 


5 where 


VarjX k ) 


(N 0 T) 



^(n,) 


, ^(n s ) + ^(n s ) +* 3 (n s > 


( 29 ) 


2 
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f exp(— 2r;,) f T)s 1' 

Cov„(X»,X t+ i) = (AfoT) [ ^ + (^T + 8, 

xerf 2 (>/r/7) + cxp(— ijj)erf(y*i7) (^J~~ + gyj 


+ ^.[f 3 (i,,) + g 3 (i),) + iM’h) - Ci(').) - M*) 

--Ft(t),)-GiM-2Gi{n,)-UhM]\ ■ ( 3 °) 

Cov 0 (X t ,X t+J ) = (pp 1 ) + H*M 

-2G H (r),)-2// s ('J.)] for J ^ 2 ( 31 ) 

The functions F if Gj and in Eqa. (28) through 
(31) are defined in Appendix B and plotted in Fig. B.l 
versus n s . Setting q s = 0 in Eqs. (20) through (31) 
yields the out-of-lock statistics in the no signal 
5 case. The no signal mean is zero but the variance is 
given by 


*«».*> ■ h aUH, (£)] 


(32) 
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(3). Square-Law Lock Detector with Nonoverlapping 
Intervals (SQNOD) 

The SQNOD detector is shown in Fig. 2b. For the 
input of Eq. (1), The inphase and quadraphase integra- 
5 tor outputs are given by 


Ik 


,(t+})T+r 

J(k + \)T+r 


r(t)dt 


f 


= < 


d k *£- + Ni(k) 
d k A{^--r) 




■+dk+lA{T-j) + Wl( fc ) 


0 < t < ? 


(33) 




r < 


T 

2 


and 


Qk = 


,(k+\)T+r 

J(k-\)T+r 


r(l)dl 


( 


= < 




dk-\A(I~-T) 

+d k A + r) + ^(k) 


0<r<£ (34) 
Z < r < Z 1 

4—^2 


The noises Nj(Jc) and N 2 (k) are given by Eqs. (15) and 
(16) after changing the integration limits to those rn 
Eqs. (33) and (34). As a result, they are independent 
Gaussian random variables with zero mean and variance 
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o z n =(N a T)/ 4. The sample X k is the difference of the 

squares (i.e., X k =ll~Ql) . From Appendix D, the in-lock 
mean and variance of Y are given by 

n sN a T (35) 


and 


nIt 2 \ 

( T|s + 24f] s +16 ^ 

M ) 

l 64 J 


( 36 ) 


5 


For the case of false lock with signal present and r 
an unknown constant over M symbols, one obtains 


Vy 0 = 0 .° 


( 37 ) 


and 

N Z 0 T 2 
M 2 j 



X 


M 


11‘ + 25ti 


60 




+M(M- 1) 


_ 2 

h s 

120 


( 38 ) 


When there is no signal present the out-of-lock mean 
is zero and the variance is given by setting r) B = 0 in 
10 Eq. (38). Consequently, the no signal out-of-lock 
variance is 
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( 39 ) 


(4). Absolute-Value Lock Detector with Honoverlapping 
Intervals (AVNOD) 

This detector shown in FIG. 3b is the same as the 
SQNOD detector shown in FIG. 3a with the squaring 
5 operations replaced by absolute value operations. 
Hence, Eqs. (32) and (33) for I k and Q k are valid but 
now X k - lljJ - |0 k | . From Appendix E, the in-lock 

statistics for Y are given by 



(40) 


and 



10 The out-of-lock mean and variance when there is a 
signal present and r an unknown constant over MT 
seconds is given by 
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»Y 0 = °-° 


(42) 


and 


o l = A [MVar a (X k ) +M(M-1) Cov o (X k ,X Kn ) ] (43) 

° m 2 


Var 0 (A%) = (No T) 


12 + 2 


~exp(-t) J ) 

4rr 


-erf 2 



1 + 3j| 


8 


0 - 


exp 



x erf 





where the function Z is defined in Appendix E and 
plotted in Fig. B.2. 

For the out-of-lock case with no signal, the mean 
5 is zero and the variance is obtained by setting q g — 
0 in Eqs. (43) through (45). Hence, the out-of-lock 
variance is given by 
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°r 0 (TU=o) 


Wl 1 _ 3 \ 
M \ 2 4n / 


(46) 


5 


(5). Signal-Power Estimator Lock Detector (SPED) 

In the detector shown in FIG. 4, denote the inte- 
grations over the first half of the assumed symbol 
interval as I k and the second half as Q k . Then, the I k 
and Q k samples are given by 


*k ~ d k 


AT 

2 




(47) 


and 


Ok = ('f _T ) +d ^ i^+N 2 (k) 


(48) 


and X k = I k Q k ' From Appendix F, the in-lock mean and 
variance of Y are 


Hr 


4 


(49) 


and 


4 = 


( 50 ) 
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The out-of-lock case with signal present has mean and 
variance when t is constant over M symbols is given by 


n sNqT 

8 


( 51 ) 


and 


<4 


NiT 2 \ 


\M 


Ar 


5ris+20g s +12 

192 




T\l 

192 


( 52 ) 


As before, the out-of-lock variance in the no signal 
5 case is given by 


< 4 (^= 0 ) 



( 53 ) 


Discussion and Simulation Results 

Digital simulation was used to verify the forego- 
ing analysis. In the out-of-lock state for a long- 
time constant B L = 1 / [MT) , the symbol timing error t 
10 is modeled as constant over a decision interval (MT 
seconds) but independent and uniformly distributed 
over the collection of all decision intervals. The 
timing error in the in-lock state is modeled as being 
zero. Although the special case of r constant over M 
15 symbols was analyzed for performance comparison pur- 
poses, it is not advisable to operate a practical 
system under these conditions due to unacceptable 
false alarm rates. This case has higher than usual 
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false alarm rates because the decision statistic for 
small values of r is not significantly different from 
the statistic for t = 0. As a result, the out-of-lock 
states corresponding to small values of r are fre- 
5 quently declared to be in-lock because they are mista- 
ken for the case when r = 0. This problem can be 
ameliorated by lengthening the observation time 
relative to the time constant of z (i.e., shortening 
the time constant of r). In practice, it is recom- 
10 mended that the observation time be at least ten times 
longer than the time constant of t. As noted above, 
the out-of-lock density function for Y in this case is 
not Gaussian. Consider the decision statistic Y when 
the loop is out of lock and r is constant over M sym- 
15 bols. In general, it can be written as 

x k = s Jt( r ) +n Jc +s Jt( r ) il k ( 54 ) 

where, in all five detectors, the signal term s k is 
random and uniformly distributed because r is a 
uniformly distributed random variable. The density of 
20 the noise n k depends on the detector being imple- 
mented. Summing M samples of X k , where X k is at the 
symbol rate in all cases, yields the decision variable 
Y. Since r is constant over the sum, at high SNR 
(i.e., for strong signal levels) the density function 
25 of Y approaches a uniform distribution as shown in 
FIG. 5a. However, at low SNR the noise term dominates 
and the density of Y is Gaussian due to the central 
limit theorem as shown in FIG. 5b. The density in 
FIGs. 5a and 5b was obtained via numerical integration 
30 as well as simulation. Both methods are seen to agree 
very well. The numerical method computed the density 
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function of Y , f (y) , by averaging over r the condi- 
tional pdf /(y/r). The latter is Gaussian with mean 
and variance both a function of z. The simulation 
method computed the histogram of Y and then set 
5 / (y) = [P(y-A<;yiy+A) ] /A , where A is the size of a 

histogram bin. The histograms were generated using 
1,00,000 million symbols which corresponds to 10,000 
decisions (F's), since there are 100 symbols/decision. 

FIG. 6 compares the probability of detection 
10 performance of all five detectors for M=100 and P j a = 
0.25. Note that the overlapping detectors SQOD and 
AVOD, which are identical except for the squaring and 
absolute value operations, have nearly identical 
performance. As expected, the AVOD is slightly better 
15 at high SNR, whereas the SQOD is slightly better at 
low SNR. The non-overlapping detectors SQNOD and 
AVNOD also have nearly equal performance. Once again, 
the absolute value operation yield better results at 
higher SNRs . The signal level estimator (SPED) is 
20 better than the non-overlapping detectors but worse 
than the overlapping detectors. The probability of 
detection results in FIG. 6 change when Pj a or M 
change. For example, increasing the observation 
interval increases the detection probability because 

25 it increases the detector SNR = (py/cry) * Accepting a 

higher false alarm rate increases the probability of 
detection because it lower the threshold 6. In 
generating these curves, 50,000 symbols were simulated 
for each value of SNR. Since there are 100 sym- 
30 bols /decision, the detection probability for a given 
SNR is based on 500 decisions. 
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In the out-of-lock state for a short-time con- 
stant ( B l - l/T) , the symbol timing error z is modeled 
as uniformly distributed random variable that changes 
independently from symbol- to-symbol . For this case, 
5 the probabilities of detection for all five detectors 
are computed by simulation for M = 100, Pj a = 10~ 2 and 
the threshold 6 is set according to Eq. (11). The 
false alarm rate was verified by simulation. The 
results are plotted in FIG. 7 versus symbol energy-to- 
10 -noise ratio q 3 . In these computer simulations, the 
detection probability for a given SNR is based on 
40,000 decisions. 

The results show that, the AVOD performs slightly 
better than SQOD at high SNR, whereas they seem to 
15 perform identically at low SNR. The nonoverlapping 
detectors SQNOD and AVNOD also have nearly equal 
performance at low SNR, but AVNOD performs about 1 dB 
better for values of symbol SNR higher than -4 dB. As 
far as the SPED, it performs about 2 dB worse than the 
20 overlapping detectors and 3 dB better than the other 
two nonoverlapping detectors. Also by simulation, the 
false-alarm rate that was used in setting the thresh- 
old was verified. 

In the situation of no signal, i.e., when there 
25 is no signal present as distinguished from the case 
when there is a signal and z - 0, the out-of-lock 
statistic is Gaussian with zero mean and the in-lock 
statistic is Gaussian with non-zero mean. Probability 
of detection results are compared in FIG. 8. Inter- 
30 estingly, the performance of the overlapping and 
nonoverlapping detectors are grouped together, but the 
signal level detector (SPED) now has the best perfor- 
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mance. The interdependence between P d , Pj a and M is 
the same as in the other two cases. 

Conclusion 

The performance of two prior-art symbol lock 
5 detectors shown in FIGs. 2a and 2b have been compared 
with three new ones shown in FIGs. 3a, 3b and 4. They 
are the square-law detector with overlapping (SQOD) 
compared with nonoverlapping (SQNOD) integrators and 
others, the absolute value detectors with overlapping 
10 (AVOD) compared with nonoverlapping (AVNOD) integra- 
tors and others, and the signal power estimator 
detector (SPED) compared with all others. The analy- 
sis considered various scenarios when the observation 
interval is much larger or equal to the symbol syn- 
15 chronizer loop bandwidth, which has not been consid- 
ered in previous analyses. Also, the case of thresh- 
old setting in the absence of signal was considered. 
The analysis has shown that the square-law detector 
with overlapping integrators (SQOD) outperforms all 
20 others when the threshold is set in the presence of a 
signal, independent of the relationship between loop 
bandwidth and observation period. The square-law 
detector and absolute-value detector with overlapping 
integrators outperformed corresponding detectors with 
25 nonoverlapping integrators, but implementation of the 
SQOD and AVOD require two separate integrators since 
both integrators must operate at the same time due to 
overlapping, whereas the SQNOD and AVNOD symbol lock 
detectors may use a single interval since the separate 
30 integration operations required do not overlap in 
time. On the other hand, the signal-power estimator 
detector (SPED) outperforms all others when the 
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threshold 6 is set in the presence of noise only, and 
it requires only a single integrator for implementa- 
tion. 
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SYMBOL LOCK DETECTION IMPLEMENTED 
WITH NONOVERLAPPING INTEGRATION INTERVALS 


ABSTRACT OF THE DISCLOSURE 

The performance of five symbol lock detectors are 
5 compared. They are the square-law detector with 
overlapping (SQOD) and non-overlapping (SQNOD) inte- 
grators, the absolute value detectors with overlapping 
and non-overlapping (AVNOD) integrators and the signal 
power estimator detector (SPED). The analysis consid- 
10 ers various scenarios when the observation interval is 
much larger or equal to the symbol synchronizer loop 
bandwidth, which has not been considered in previous 
analyses. Also, the case of threshold setting in the 
absence of signal is considered. It is shown that the 
15 SQOD outperforms all others when the threshold is set 
in the presence of signal, independent of the rela- 
tionship between loop bandwidth and observation 
period. On the other hand, the SPED outperforms all 
others when the threshold is set in the presence of 
noise only. 


20 
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Appendix A 

Derivation of the Mean and Variance of the SQOD 

The inphase and quadraphase integrator outputs are given by Eqs. (13) and (14), respectively. The output of the 
lock detector A* = I k ' - Q*. Consequently, 




(A-l) 


where ft k = £{X k ). 


Var(*0 = £{[/*’- Q*T) - - QZ) 

^eVt + Qt-utQD-vl 


(A-2) 


The covariance of X k with X k+i is 


Cov(* 4 ,* i+> ) = £ {1W +1 ) + £{QlQl +j ) - £{llQl H ) - WUjQl) - Hk „ l+i 


(A-3) 


When the loop is in lock, Eqs. (A-l) through (A-3) are evaluated with r set to zero in Eqs. (A-10) through (A-18) 
Hence, the in-lock moments of X* are given by ' b V '* 


A 

Var(Afi) = — — -f 8A 2 T 2 a 2 + 12 a* 


(A-4) 


and 


Cov(X k ,X k + 1 ) = -2A 2 T 2 <r 2 n -2* i n 


(A-5) 


and for j > 2, this can be shown to be 


Cov(X k ,X k+j )=0 


(A-6) 


th ™ h wi h ! H ^ > ° P ° l ut ; of - ,ock - l is Celled as a uniform random variable. Using this model for r in Eqs. (A-10) 

ler t. ni r S | i ,tUt,,,g th , e reSUlU int ° EqS - ( A ' ] ) lhrou fi'* (A-3) give the out-of-lock moments of X k . Namely 
(where the additional subscript o denotes out-of-lock), 


v».« - Ob = ^ + 


~ , v v v 23A 4 7 V * 23 , 

Cov e (A t , A t+1 ) = -* 32 q- - j 2 A T a l ~ 


(A-7) 

(A-8) 


and for j > 2, this can be shown to be 
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Cov 9 (X k ,X k+i ) = (A-9) 

The following equations were used to compute the variance of X k and covariance of X k with A’* + 1 : 


£{ll) = A\T 2 - 2 T£{r) + 2£ {r 2 }) + 2<r 2 (A-10) 

, (All) 

{ A 2 (Z£-€ 3 {1) -4T£ 3 {r) + 2£ 2 {r 2 )) + 2<r 2 £ 2 {l) |<r<T 


£{li) = A'(T* - 4T 3 £{r) + 12T 2 £{r 2 } - 16TC{r 3 } + 8£ {r 4 }) + A 2 <r 2 (12T 2 - 24 T£{r} + 24 £{r 2 )) + 12a< 


+ 8 ^i{r < }) 

+A 2 <T 2 (6r 2 f l {l} + 24A{r 2 }) + 12 < r^ l {l) 0<r< £ 

/t<(i^£ 2 {l} - 32T 3 £ 2 {r} +48T 2 £ 2 {r 2 } - 32T^{r 3 ) + 8£,{r 4 }) 

. +^(30T 2 f a {l}-48Tf 2 {r} + 24r 2 {r 2 }) + 12a^ 2 {l} Z < r < T 


(A- 12) 


(A- 13) 


£{ilQi) = { 


A A (^-£ t { 1} + 2T 2 £, {r 2 } - 8 T£ x {r 3 } + 8£, {r<}) 
+/l 2 a 2 (5r 2 f,{l}-8Tf,{r)+16A{r 2 }) + 6«r^ I {l} 
A‘*(^ 2 {l}-9T 3 f 2 {r}+ I5r 2 f 2 (r 2 }- 12Tf ? {r 3 ) + 4f 2 {r < }) 

+A 2 ff '-(7T i £ 2 {l}-6T£ 1 {r) + 4£ 3 {T i )) + 6a i n £ 2 {\) 


0 < r < I 
T < r < T 


(A- 14) 


£{I\ Uhl) = - 4 T 3 £(t) + 8T 2 £{t 2 ) - 8T£ {r 3 } + 4£{r<}) + /t 2 <r 2 (4T 2 


-8T£:{r)+8r{r 2 }) + 4^ 


(A- 1 5) 


f ^(^ i {l)+7T 2 £ i (T 2 )+4£ l {r A )) 


^{QlQ2+,} = { 


+A 2 <r 2 (2T 2 A{l)+8f,{r 2 )) + 4^f l {l) 

- 20T 3 £ 2 {r} + 26T 2 £ 2 {r 2 } - 16r£ 2 (r 3 } + 4£ 2 {r'}) 


0 < t< \ 


(A- 16) 


l + /» 2 ^(IO'f 2 £ 2 {l) - I6T£ 2 {r) + 8£ 2 (r 2 )) + 4 <t*£ 2 { 1 } $ < r < T 
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' ^(^,{1} - T 3 r,{r} 4 ZT'rAr'*) - 4Tf,{r 3 } 4 ^,{r‘*}) 
4 A j <t 2 (3T j £, {1 } - 4T£i{t) + 8£, {r 1 }) + 4a^i {1} 

' ^(5Tlf 2 {l) - 9 T 3 £ 2 {t) + 15 T'&Ar 7 ) - 12 TC 2 {r 3 ) 4 4£ 2 {^1) 
+A'ff*(7T 7 ei{\) - UTC^r) 4 8f 2 {r l })4 4a^ 2 {l} 


0 < T < 7 
j<r<T 


A\^£i{l) -T^ir) 4 3T J f,{r 2 ) -4T£:,{r 3 } 4 4A{^}) 


£{iUtQl) = { 


+A*<rl(bT 7 £i{\) - 2TA{r} 4 4^,{r 2 }) 4 6 <t^i{ 1) 
yt«(Sf £,{1} - 22lif 2 {r) 4 ^{r*} - 15Tf 2 {r 3 ) 4 5£ 2 {t<}) 
4/lV 2 (13T 2 f 2 {l} -24T5 2 {r} 4 16 £ 2 {t 2 }) 4 6<r<£ 2 {l} 


0 < r < y 
y < r < T 


where, in the above equations 


£{/(t)} - / /( r )p( r ) <fr 
Jo 

£j{/(r)} = / /(r)p(r)dr 
Jo 

£ 2 {/(r)} = [ T nr)p(r)dr 
•'T 


(A- 17) 

(A- 1 8) 

(A- 19) 
(A-20) 
(A-21) 


where p(r) is the probabilty density function of the variable r. 



- 39 ^ 


Appendix B 

Derivation of the Mean and Variance of the AVOD 


Note that the calculations in this appendix incorporate the results of Appendix C. The inphase and quadraphase 
outputs are given, respectively, by Eqs. (13) and (14). The lock detector output X * = |/*| — |Q*|. Let /j* = £ { | /jt | — |Qt|}- 
Then, 


Var(A*) = f{[|/ t |-|C? t |] 1 }-^ 


= £{ll + Ql- 


(B-l) 


and 


Cov(* t ,* fc+ i) = £{[\h I - |Q* I - PkWh+A ~ IQk+il -w+i]} 

s^|/*/* +i |} + f{|Q*Q* +j |} -f{|/*Q*«|} -f{|/*+jQ*l) -/•*/**+> ( D -2) 


The following equations were used to compute 


the variance of X* and the covariance of X* with Xk+\' 


( \{£ X {| AT - 2 At + Af |) - £ x {|2 At + Af|)) 0 < r < \ 

£{X k ) = (B-3) 

{ i[£ 2 (|AT-2Ar+N|}-£ 2 {|2AT-2Ar+VV|}] £ <r<T 


£{\hQk \ } = 


r i(A {|(AT + N x (k) + 7V 2 (*)][AT + N 2 (k) + N x (k + 1)]|) 

+£ x {| [AT - 2 At + N x (k) + Af 2 (fr)](-2Ar + Af 2 (Jt) + N x (k + 1)]|}] 
\{£,{\ [AT + Ni(k) + N,(k))[AT+N»(k) + JV x (k + 1)J|) 

+£ 3 {UAT+ Ni(k) + H 2 (k)][2AT-2AT + N 2 {k) + N x (k + 1))|) 

+/ 2 {| [AT - 2At + ATi(Jt) + Nt{k)][2AT-2Ar + N 2 (k) + N x (k + 1)]|) 
s +£ 2 {\[AT - 2Ar + N\(k) + Af,(*)pT + l V 2 (i) + N x (k + 1)J|}] 


0 < r < y 


y < r < T 


(B-4) 


£{|/*/*+i|) = \[£ 7 {\AT+ N x (k) + N 2 (k)\) 

+ £{\ [AT - 2 At + N,(Jk) + N 2 (k)][AT - 2 At + N x (k + 1) + JV 2 (Jk + 1))|) 

+ 2£(\AT + A/,(Jt) + W 2 (Jfc)|}£{|AT- 2Ar + N x (k + 1) + N 2 (Jk + 1)|}] (B-5) 
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' N x + AT 2 |} + f, {|[2/lr + Nt(k + 1) + W 2 (*)] 

x [2At -f N\(k + 2) -f N 2 (k + 1))|) 

+2£{\AT + TV, + N 2 \)S x {\2At + N x + at 2 |}] 

\[£^{\AT + N { + 7V 2 |) + €i{\[2AT - 2 At + AT, (it + 1) + N 2 (k)] 
x [2AT - 2At + N t (k + 2) + N 2 (k + 1)J|} 

+2£{\AT + N X + N 2 \)£ 2 {\2AT - 2 At + AT, + JV 2 |}] 


0 < r < i 


(B-6) 


j<r<T 


£{\hQk+i\) 


' \\^{\AT-kN x -f 7V 2 |} + £ l {|[2Ar + JV,(Jb + l) + yv 2 (jl-)j 
x[AT - 2 At + N x (k + 2) + N 2 (k + 1))|} 

+2£{\AT+ N> + N 2 \)£ 1 {\2At + N x + JV 2 |}) 0 < r < £ 

\VH\AT + N x + W 2 |} + £ 2 {\[AT - 2At + N x (k + 1) + Af 2 (Jt)] 
x [2 AT - 2 At + N x (Jt + 2) + N 2 (k + 1)]|) 

. +2£{\AT+ N x + N 2 \)£{\AT-2Ar+ N t +N 2 \)] j<r<T 


(B-7) 


f{|7 t+ ,Q t |} = { 


f ?(^{|[4r+^,(t + l) + ^ 2 (Jt+ l)pr+ A r i(Jt + l)+yV 2 (i)]|) 

+A {|[/tT - 2At + Ni(k + 1) + W 2 (* + 1)][/17’+ Af,(* + 1) + // 2 (*)]|} 
+£i{\[AT - 2Ar + N x (k + 1) + N 2 (k + \)\[2Ar + N x (k + 1) + AT 2 (jt))|} 

+£, {| [AT + N x (k + 1) + N 2 (k + 1)]|2 At + N x (k + 1) + W a (Jb)J|}] 
i[£?{\\AT-2At + N x {k+ \) + N 2 (k + l)][2A7 , -2Ar + Ar,(Jb + 1) + // 2 (ib)]|} 
+£ 2 (\[AT + N x (k+ 1) + N 2 (k + I)p7’+W,(H 1) + W 2 (Jb)]|}] 


0 < T < y 


(D-8) 


y < r < T 


where £ x and £ 2 are defined in Appendix A. The following functions were defined to obtain the results in Subsection 
III. A of the main text. 


= £{|(1 + cnj + cn 2 )(l +cn 2 + cn 3 )|} 

Fi(*h) = £{|(1 + cnj + cn 2 )(cn 2 + cn 3 )|) 

Fs(n>) - £{|(1 - 2u + cn! + cn 2 )(l - 2u + cn 3 + cn 4 )|} 
= £i{|(] — 2u + cnj + cn 2 )(2u + cn 2 + cn 3 )|} 
G 2 (*u) = £\ { |( 1 - 2u -f cni + cn 2 )(l + cn 2 + cn 3 )|} 
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GsM = A {|(2u + cm + cm)(2u + cn 3 + cn 4 )|) 

G 4 (rj,) = £1 { |(2u + cm + cn 2 )(l - 2u + cn 3 + cn 4 )|} 

// l (rj J ) = 5i{|(l + cm + cn 2 )(2-2u + cnj + cn 3 )|} 

tf 2 (rj,) = £{|(1 - 2u + cm + cnj)(2 - 2u + cn 2 + cn 3 )|) 

HM = £{|(l - 2u + cm + cm)(i + cm + cn 3 )|} 

H 4 (q,) = S {|(2 - 2u + cm + cm)(2 - 2u + cn 3 + cn 4 )|} 

= £{](1 - 2u + cm + crn)(2 - 2u + cn 3 + cn 4 )|} (B-9) 

where c = l/(2.i/rj7), the n,’s are normal independent random variables with zero mean and unit variance, and u is a 
uniform random variable in the range [0, 1]. In Fig. B-l, one plots these functions versus n . These functions have been 
computed as follows: In the F functions, the expectation with respect to u is earned over the entire region (0, lj, while 
in G and H functions, the expectation is carried over [0, 1/2] and (1/2, 1], respectively. 





Appendix C 

Derivation of the Mean for the Random Variable \n + bt + c\ 


Let n be a normal random variable with zero mean and variance a 2 , r a uniform random variable over (0, T), and c 
any constant. Then, 

✓ . 2 \ 

(C-l) 


fi|n + = vh £ |n + c, “ p (■£) 


The above integral can be easily evaluated by breaking it into two integrals over the two regions (- 00 , -c) and ( c, 00 ) 
to obtain 


£{\ 


” + .ii=\/?' xp (^) +ce,f fe) 


(0-2) 


For a fixed r, one can write 


€{\n + br + c\) = /^exp + (br + c)erf 

Unconditioning over r yields 

€i{\n + bT + c\} = j;J* [\/v exp ("^ 2 ^" ) + (6r + c)erf ( 7 ^ 1 ) 

which, after integration by parts, leads to 
£»{|n + br + c|} = 

^ ^ V ~ ] \ V * 

*S[(¥-)-(=^)--( iS)] 

By applying the above expression, one gets 

S\ {\AT - 2At +N,+ Nj\) = |^=erf( Ju) + ~tr((^nT} + exp(“'?-) 

Also, by simple manipulation, it can be shown that 

£ 1 {\AT-2At + N 1 + N 2 \) = £ 1 {\AT-2At + Ni + N 7 \) 

= £ { {\2At + Ni + N 7 \) 


dr 


<r 2 


erf ( c \ 

1 

4- - - 

f- + C Verf( 

'•F + A 

— c J erf ( l-') 

2bT 

rwJ 

_erf VV2^J 

2bT 

\2 + ) \ 

^ Vtt ) 



(C-3) 


(04) 


(C-5) 


(C-6) 


= £j{\2AT - 2Ar + AT, + A/ 2 |} 


(0-7) 



Appendix D 

Derivation of the Mean and Variance of the SQNOD 


The samples h and Q k are given by Eqs. (33) and (3d). The moments of A'* are given by Eqs. (A-l) through (A-3). 
Using Eqs. (D-7) through (D-10) with r = 0 and Eqs. (A-2) and (A-3) yields the in-lock variance. Namely, 

v.,(X,) = ^p + |/’r>.; + <»l (D>) 


The covariance can be shown to be, for j > 1, 


Cov(A'*,Af t+i ) = 0 (D-2) 

Similarly, with r modelled as uniform over (0, Tj 2) in Eqs. (D-7) through (D-10), the out-of-lock moments are found. 
Hence, 


Var.(A' t ) 


A*T* 

CO 


+ Ia't’vI + u 


4 

n 


The out-of-lock covariance can be shown to be, for j > 1, 


Cov.(X k ,X k+j ) 


A 4 T< 

120 


(D-3) 


(D-4) 


The following equations are used in computing the variance of X k : 

£{!}) = 

U 2 (^ 

» 2 (£Ml) + 2£.{r 2 )) + o^,{l) 0<r<£ 


0 < T < j 


+ i) 

A 2 {^-£ 2 {1) -2T£ 2 {t) + 2£ 2 {t 2 )) + * 2 „£ 2 {1) £ <r< \ 


e[Ql) = 




Z < r < Z 

4-^5 


, (4^1+|^r 2 < r 2 + 3^K l {l} 

£{lt) = < A\yg-£i{ 1} -4T^^ a <r) + UT 2 £ 2 {t 2 ) - 16T£ 2 {t 3 } + 6£ 2 {t a )) 
, +A 2 (rlCiT 2 £ 2 { 1) - \ 1 T£ 2 {t) + U£ 2 {t 2 )) + ic*£ 2 {l) 


0 < r < 


Z < T < Z 

4 — T ^ 2 


£{Qi) = 


j A\^£i (1) + 8£,{r<}) + A 2 a 2 n {\T 2 £ k (\) + 12£,{r 2 )) + 2a*£ t {1} 0 < r < £ 

l(^ + ^ 2 TV n } + 30^1) 


Z < r < T- 

4 - T ^ 2 


(D-5) 


(D-G) 


(D-7) 


(D-8) 



ijch 
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Appendix E 

Derivation of the Mean and Variance of the AVNOD 


By following Ihe same procedure as in Appendix B and by using 


£{|/*n = 


$(£ 2 {|4f + N\) + £i{\AT - 2 At + AM}] 


0<r<^ 
T — r < 7 


(E-l) 


£{\Qk\) 


^,tl^ + Af|}+A{|2A r + Ar|}] o <r<^ 


(E-2) 


and 


Var(|/ t |-|Q*|) = 


' 2dlTifi{l} + 2A J f,{r J } + 2 < r*A{l} 

-|A J {|^ + ^1}-K?{|2^ + Af|} 

-|f{|^ + Af|K,{|2Ar + /V|} 
l££-£l{l) -2A*T£- l {r) +2A 7 £i{r 7 ) +2al£ 7 [l) 

-f^{|^ + Ar|)-Klil AT - 2Ar + // li 

k -\£{\Af + N\)£ 2 {\AT-2Ar + N\) 


0 < r < j 


k< 


T < 


T 

T 


(E-3) 


one obtains Eqs. (40) through (45) after using the results of Appendix C and lengthy manipulations. The function Z in 
Eq. (45) is defined as 


Z(rj,) = £i{|(u + cni)(u + cn 2 )|} 


(E-4) 


where c = l/(2^ri7), the n/s are normal independent random variables with zero mean and unit variance and u is a 
uniform random variable in the range (0, 1). In Z, the expectation with respect to u is over the range [ , / 1 ai ' 

function is plotted in Fig. B-2 versus the symbol SNR (tyj). 


Appendix F 

Derivation of the Mean and Variance of the SPED 

raightforward to show that" 6 g Ve " ^ ^ ^ The out P ul of the lock detector is 


riy 1 — A^T* A*T 
£{X k )= £{r) 


and 


x * — hQt . It is 


(F-l) 


£{X 7 k ) = a* (— - Tl£tl , T*£{r') \ , T 1 

VI® * 2 ) + ff » + A < (- 5 - ~ TC{r) + 2£{r i ) j 


Cov(A' t , A' i+J ) s J ~ Var ( r ) r ~ “(0. 7) a nd > > 1 


r = 0 


(F-2) 


(F-3) 


r be unifo/m iver [O^T/^ ^ T ° ,n EqS ( F * 1 ) trough (F-3). Equations (51) and (52) follow by letting 
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